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Abstract 

In a typical three-dimensional metal, the low-energy excitations are found on a two-dimensional 
closed Fermi surface in momentum space. Topological semimetals, by contrast, can support one¬ 
dimensional Fermi lines or zero-dimensional Fermi-Weyl points, at locations in momentum space, 
where the valence and conduction bands touch. While the degeneracy points in Weyl semimet¬ 
als are protected against any perturbation that preserves translational symmetry, nodal lines are 
topologically protected by symmetries such as mirror reflection. Here we report on the existence 
of topological nodal-line states in the non-centrosymmetric compound single-crystalline PbTaSe 2 
with strong spin-orbit coupling. Remarkably, the spin-orbit nodal lines in PbTaSe 2 are not only 
protected by the reflection symmetry but also characterized by an integer topological invariant as 
we show here. Our detailed angle-resolved photoemission (ARPES) measurements, first-principles 
simulations and theoretical analysis illustrate the physical mechanism underlying the formation 
of the topological nodal-line states (TNLS, in contrast to Weyl Semimetals) and associated sur¬ 
face states for the first time. Therefore, this work paves the way towards exploring the exotic 
properties of the topological nodal-line fermions in condensed matter systems and, potentially, the 
rich physics arising from the interplay between the topological nodal-line states and the emergent 
superconductivity we report in the single crystal phase phase of this compound. 
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The discovery of the time-reversal invariant topological insulator has stimulated enor¬ 
mous research interests in novel topological states protected by symmetries [i[J^ . One of 
the key properties of topological materials is the existence of symmetry-protected metallic 
edge or surface modes in bulk-insulating ground states which is due to a topologically non¬ 


trivial ordering o: 
Weyl semimetals 


Dulk wave functions. Recently, due to the experimental observations of 
4l-lll| , the research interest in topological phenomena in condensed matter 
systems has partially shifted to from insulators to semimetals and metals. A Weyl semimetal 
is a topological state of matter whose low-energy bulk electrons are linearly dispersing Weyl 
fermions. The two-fold degenerate Weyl nodes, carrying non-zero chiral charge, are con¬ 
nected on the boundary by Fermi arc surface states, which are predicted to exhibit unusual 


transport behaviors 12Hl5|. In contrast to Weyl semimetals whose bulk Fermi surface has 


dimension zero, topological nodal-line semimetals have extended band touching along one¬ 
dimensional curves in k space, presenting a significant expansion of topological materials 
beyond topological insulators and Weyl semimetals, and new opportunities to explore exotic 
topological nodal physics. Line-like touchings of a conduction and valence band require ex¬ 
tra symmetries besides translation, such as mirror reflection, to be topologically protected. 
Kinematically this protection involves a finite fraction of Brillouin zone. For this reason, thi s 
leads potentially to many anomalies in their electromagnetic and transport response 


16 


181 . 


Similar to the case of Weyl nodes, one can define an i nteg er topological invariant for the 
line node along which two nondegenerate bands touch 1^. Despite the many theoretical 
discussions of nodal-line semimetals, a material realization of topological nodal-line fermions 
has been lacking for many years just like Weyl semimetals. Here we report the existence of 
topological nodal-line phase in the normal states of the superconducting compound PbTaSe 2 
for the first time. 

The crystal lattice of PbTaSe 2 lacks space inversion symmetry, which lifts the spin- 
degeneracy of its electronic bulk bands. Our angle-resolved photoemission (ARPES) mea¬ 
surements together with DFT calculations show that the conduction band originated from 
Pb-6p orbitals and the valence band from Ta-Sd orbitals cross each other, forming three 
nodal-line states close to the Fermi energy. The nodal lines are protected by a reflection 
symmetry of the space group as we show. The topological-nodal-line state in PbTaSe 2 be¬ 


longs to the symmetry class A+R {p = 2) of symmetry-protected semimetals Q . We also 
demonstrate through effective Hamiltonian modeling and DFT simulations that the nodal 
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lines are accompanied by unique surface bands. These topological surface states are due 
to the TT Berry phase agglomerated around the nodal line in analogy to the states on the 
graphene-like zigzag edge. Our experimental and theoretical results establish a material 
realization of topological nodal-line fermions in the superconducting compound PbTaSe 2 . 
The hndings reported here are of substantial importance because the door for exploring the 
exotic properties of nodal-line states in condensed matter and, potentially, the rich physics 
arising from the interplay of nodal lines and emergent superconductivity has been opened 
by the identihcation of this compound. 

Our PbTaSe 2 single crystals were prepared by the chemical vapor transport (CVT) 
method, see Fig. la. The samples were of high structural quality, which was conhrmed 
by our X-ray diffraction (XRD) and scanning tunneling microscopy (STM) measurements. 
The XRD peaks shown in Fig. lb are consistent with the space group of PbTaSe 2 , PQm2 
(187), and, therefore demonstrating the lack of inversion symmetry of our PbTaSe 2 single 
crystals. This is crucially important for lifting the spin degeneracy, a necessary condition 
for the formation of topological nodal lines. To further check the chemical composition of 
our samples, we performed a photoemission core level scan. Clear Pb-Sd, Ta-4/ and Se-3d 
core level peaks were observed in the photoemission spectrum, which conhrms the correct 
chemical composition in our PbTaSe 2 single crysta samples, shown in Fig. Ic. To verify the 
superconducting property of our samples, a transport measurement was carried out. The 
measured resistivity curve. Fig. Id, shows a clear superconducting transition temperature 
at 3.8 K, consistent with the value reported in 19(]. Figures le and If show STM images 
of the cleaved (001) surface. The topography image clearly reveals a hexagonal lattice with 
few defects, demonstrating the high quality of our samples. Furthermore, no surface recon¬ 
struction was observed on the cleaved surface. The high-resolution STM topography yields 
a lattice constant of 3.2 A. 

PbTaSe 2 crystalizes in a hexagonal lattice system in which the unit cell consists of one 
Pb, one Ta and two Se atoms and each atom resides on a hexagonal flat layer, shown in 
Fig. 2a. The stacking sequence of these atomic planes within the unit cell is Pb-Se-Ta-Se: 
A-A-B-A (A, B and C, here, refer to the three high-symmetry spots on a hexagonal lattice). 
The lattice can also be viewed as a Pb-layer intercalating two adjacent TaSe 2 layers with 
Pb atoms sitting above Se atoms. The Pb intercalation suppresses the softening of phonon 
modes associated with the charge density wave (CDW) in TaSe 2 and stabilizes the hexagonal 
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lattice on the surface 


20| . This particular stacking does not preserve the space inversion 


symmetry, but the lattice is reflection-symmetric with respect to the Ta atomic plane. In 
other words, the Ta atomic planes are a mirror plane of the crystal lattice under the mirror 
operation that sends .2 to —z. This reflection symmetry of the lattice provides a protection 
for the topological nodal lines, as discussed later on. The bulk and (OOl)-projected surface 
Brillouin zones are shown in Fig. 2b. The A, H and L points are high symmetry points on 
the = IT plane, which is a mirror plane of the bulk Brillouin zone. Figure 2c presents 
an overview of the band structure calculation for PbTaSe 2 , which was performed by the 
DFT-GGA method. Glose to the Fermi level, two prominent features in the band structure 
are observed. One is a giant hole pocket around F, whose states are mainly derived from the 
Ta 5(i322_r2 orbitals that are oriented out of the Ta atomic plane, taking the Ta plane as the 
x — y plane. The second major contribution to the density of states at the Fermi level comes 
from the four bands that cross each other near H. The two electron-like conduction bands 
originate from Ph-6px/py orbitals and the two hole-like valence bands from Ta-bclxy/dx^-y^ 
orbitals. We note that all these orbitals are invariant under Rz. A zoom-in view of the bands 
around H without/with spin-orbit coupling (SOG) is shown in Fig. 2d and 2(e), respectively. 
Without the inclusion of SOG, the conduction and valence bands become spin-degenerate. 
The two bands belong to different representations of the space group (the representation of 
the electron-like band is A' and that of the hole-like band is A"), therefore the intersection of 
the two bands is protected by the crystalline symmetry, forming a spinless nodal ring. Once 
SOG is turned on, each band split into two spin branches with opposite spin orientations and 
mirror reflection eigenvalues as indicated in Fig. 2e. Only the crossings of branches with 
opposite mirror reflection eigenvalues remain gapless as a result of symmetry protection, 
forming a pair of nodal rings. Interestingly, SOG also gives rise to a third nodal ring on 
^2 = 0 plane. The detailed band dispersion and the rise of three nodal rings are very well 
captured by our effective k • p Hamiltonian, please see |2l|. Before proceeding to a detailed 
discussion of the nodal-line states, we will present the results of our ARPES measurement, 
verifying the overall band dispersion of Pb conduction bands and Ta valence bands obtained 
from our DFT calculation. 

Fig. 3a shows a brief overview of our ARPES band mapping and the corresponding 
numerical calculation of the PbTaSe 2 band structure is presented in Fig. 3b. The projected 
bulk bands and surface bands (as highlighted by white lines) were calculated for the Pb- 
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terminated (001) surface. The DFT band structure reproduces the ARPES spectrum very 
well. Specihcally, in the ARPES spectral cut, a band marked as SSi with high intensity 
poke the Fermi level between F and K. This is the surface state band associated with the 
Pb-terminated (001) surface 2l(]. Around K there are three concave bands whose binding 
energy at K are 0.21, 0.75, and 0.80 eV, respectively. The top and bottom bands correspond 
to the electron-like bands derived from Pb-6p orbitals. The middle band, marked as SS 2 , is 
consistent with the surface band as plotted in Fig. 3b. The two bands at M are tails of the 
two Ta-5d bands that cross the two Pb-6p bands forming the nodal rings in the vicinity of 
K. Two Ta-5d bands have to degenerate in energy at M according to the Kramers theorem. 
The ARPES measured (001) Fermi surface with the incident photon energy of 64 eV and 
the theoretical simulation are shown in Figs. 3c and 3d, respectively. At the Fermi level, our 
data shows that the Fermi surface consists of three parts: a hexagon-shaped pocket centered 
at F with smeared intensity inside, a dog-bone shaped contour centered at the M point and 
several circles surrounding the K point. Our ARPES data and calculation show agreement 
on those features. Furthermore, the hexagon centered at F and the intensity inside are the 
surface band and the bulk hole pocket at F, respectively. The dog-bone shaped contour 
corresponds to the one branch of the Ta valence band and the circles around K are from 
the other branch of the Ta valence band, the surface states and the spin-split conduction 
band derived from Pb orbitals. As the binding energy decreases we hnd that the Pb pockets 
at K shrink while the Ta pockets expand outwards, which is in good accordance with the 
characteristics of the electron-like Pb bands and hole-like Ta bands. 

TaSe 2 can be regarded as a building block of PbTaSe 2 , and therefore its electronic struc¬ 
ture can be traced from that of PbTaSe 2 . To highlight the difference between electronic 
structure of TaSe 2 and PbTaSe 2 , we mapped out the Fermi surface and band structure 
along M-K-M of the two compounds, shown in Figs. 4a and 4b. In the Fermi surface map¬ 
ping of TaSe 2 , there are one dog-bone shaped contour centered at M and only one circle 
shaped contour centered at K. Those contours are from Ta valence bands, and consistent 
with previous work 22|. By contrast, the Fermi surface of PbTaSe 2 has more ring-shaped 
contours centered at K, signifying the contribution from the Pb layers. It is easier to view 
this difference from the M-K-M cut. TaSe 2 does not show any electron-like bands at K that 
exist in PbTaSe 2 . Fig. 4c shows the ARPES mapping of the Pb and Ta bands of PbTaSe 2 
with photon energies from 54 eV to 70 eV. The middle band at K does not show any photon- 
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energy dependence, which is consistent with the surface nature of this band. However the 
other bands at K and M don’t exhibit obvious changes with different photon energies either. 
This seems to contradict to the assignment of those band as bulk bands according to our 
DFT calculations. The inconsistency can be understood by considering the fact that the Pb- 
Qpx/py and Ta-bdxy/dx^-y^ orbitals that constitute those bands are primarily conhned within 
the Pb and Ta atomic planes (which are parallel to the x-y plane), and, thus, the interlayer 
couplings (say, the coupling of one Ph-6px/py orbital with another orbital on the adjacent 
Pb plane) is largely suppressed, which results in little /photon-energy dependence. 

From the discussion before, we know that the electron-like bands from the intercalated Pb 
layers are the essential component for forming the topological-nodal-ring band structure. By 
comparing with the TaSe 2 spectrum, our ARPES established unambiguously the existence of 
the Pb bands. To further examine the topological-nodal-line states and associated surface 
states, we calculated the band structure for Pb- and Se-terminated surfaces as shown in 
Figs. 5a-d. The projected bulk band on each cut from K shows three nodal points at 
0.05 eV, 0.15 eV, and 0.03 eV above the Fermi level. The hrst two closer to K lie on the 
^2 = tt plane while the third one is on the = 0 plane. Let us refer to these three nodal 
lines as NLl, NL2 and NL3. Corresponding nodal points can be found on a cut of arbitrary 
orientation that includes K. For example, the band structure along a generic direction K—X 
is shown in Figs. 5b and 5d. Unlike the projected bulk band which is independent of surface 
termination, the dispersion of surface bands is found to be sensitive to the surface condition. 
However, in both cases we do hnd a surface band connecting to each nodal line, indicative 
of the topological nature of the bulk nodal lines. In the Pb-termination case, the surface 
bands disperse outwards with respect to K, from NLl. The surface band connecting to NL2 
grazes inwards at the edge of the lower bulk Dirac cone and merge into the bulk band. The 
surface band from NL3 disperses inwards with respect to K, consistent with the “SS 2 ” band 
in our ARPES spectrum in Fig. 3, which forms a “drumhead” surface state contour. By 
contrast, on the Se-terminated surface the surface band connecting to NL2 hrst moves into 
the bulk band gap and then fall into the bulk band region. The surface band from NLl 


disperses outwards and connects to NL3. Please refer to 


2l[ | for a detailed visualization 


of connection of the surface bands to the bulk nodal lines. To get an overall view of the 
nodal ring and surface band, we plot in Fig. 5e the isoenergy contour in the vicinity of the 
NLl nodal ring of the Pb-terminated surface and the NL2 nodal ring of the Se-terminated 



surface, as indicated by the red dashed lines in Figs. 5a-d. Indeed, gapless nodal points and 
surface states can be found at every in-plane angle departing from K. These nodal rings are 
protected against gap opening by the crystalline symmetry. Specihcally, the states in the 
two Pb branches belong to two different representations of the space group, namely S 3 and 
S 4 as shown in Fig. 5f. The same is true for the two Se branches. In particular, with respect 
to the Ta atomic plane, the two representations have opposite mirror eigenvalues under the 
reflection operation. Therefore, gap opening is forbidden at the crossing point between two 
branches of different mirror eigenvalues, which results in the nodal rings discussed in this 
work. In this sense, the nodal rings is under the protection of the reflection symmetry. If we 
shift the Pb atom slightly in the vertical direction thus breaking the reflection symmetry, all 
of the four branches are found to belong to the same S 2 representation of the reduced space 
group and, in this case, a gap opening is allowed at every crossing point of these branches 
as illustrated in Fig. 5g. A similar gap opening is also found in NL3 on = 0 plane upon 
breaking the reflection symmetry 21 ]. 

Let us briefly discuss the topological characterization of the nodal lines and the origin 
of the surface bands. The material PbTaSe 2 is time-reversal symmetric, with time-reversal 
symmetry is represented by T = cr 2 /C, where K. denotes complex conjugation and <72 is the 
second Pauli matrix acting on the electron spin. The mirror symmetry acts in spin 
space as and therefore anticommutes with T. This would place PbTaSe 2 in class AII-i?_ 
in the classihcation of [l 6 |. However, since the nodal lines are centered around momenta 
H/H' and K/K' which are not invariant under time-reversal, but pairwise map into each 
other, the time-reversal symmetry imposes no constraints on the nodal lines individually. 
The material has therefore to be classihed according to the time-reversal breaking class A-i? 
which admits a integer topological classihcation for Fermi surfaces of codimension 2, i.e., 
lines (p = 2 in [^). The nodal lines carry a topological quantum number n"*", which is given 
by the difference in the number of occupied bands with eigenvalue -|-1 inside and outside 
the line. In the case at hand, rR' = —1 for the nodal line (NL3) in the ^2 = 0 plane, while 

= -|-1 (NLl) and rR' = —1 (NL2) for the two nodal lines in the plane. We have 

also computed numerically under the DFT framework the winding number 7 = ^ dfc • A{k), 
where A{k) = is the Berry connection of the occupied Bloch bands jua.fc)- 

For a closed loop encircling any of the nodal nodal lines, we hnd that 7 = ± 7 r with the same 
sign as n'^ of that line, as shown in Fig. 5h. 
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The topological origin of the observed surface states is rather subtle. Surface states 
associated with the topological invariant via the bulk-boundary correspondence should 
only appear on surfaces that preserve R^- The (001) surface, however, breaks R^- The 
reason why we still observe surface states can be understood from the Berry phase of vr 
around the nodal line and the analogy to the edge states on the zigzag edge of graphene. 
Consider the bulk Hamiltonian on a plane in momentum space that contains both K and 
H. At low energies, each nodal ring pierces this plane twice giving rise to two Dirac cones 
in this Hamiltonian. These two cones have a Berry phase Rtt and —tt with respect to the 
orientation of the plane, precisely as in (spinless) graphene, as schematically depicted in Fig. 
5i. We know that for any termination of a graphene sample, an edge state emanates from the 
projection of each of the Dirac points in the edge Brillouin zone (except for the pathological 
case where both Dirac points project on the same spot in the edge Brillouin zone). By 
this analogy, we also expect these edge states to emanate from the surface projections of the 
nodal lines in any direction away from the K point, thereby forming a surface band. We note 
that the dispersion of this band and even whether it appears inside or outside the projection 
of the nodal line is not universal and depends on the details of the surface termination. 

Recently, some preprint theoretical works have reported on Ca 3 P 2 and CusPdN proposing 
there may be nodal-line states 23|, [2^ . We note that our work is distinct from those works 


in three aspects: (1) Both Ca 3 P 2 and CusPdN are centrosymmetric and, therefore, due to 
the coexistence of time-reversal and inversion symmetry, have four-fold degeneracy at the 
nodal ring. By contrast, the degeneracy of nodal-ring states is two in PbTaSe 2 due to the 
lack of inversion symmetry. (2) In Ca 3 P 2 and CusPdN, the nodal-line states exist only in 
the absence of SOC. In real materials SOC, however, always exists. In PbTaSe 2 , SOC is 
essential ingredient for the formation of nodal-ring states. (3) Our ARPES measurement 
established an experimental characterization of the topological nodal-line material PbTaSe 2 . 

In summary, our direct experimental observation by ARPES of the coexistence of Pb 
concave bands and Ta convex bands centered at K in the non-centrosymmetric supercon¬ 
ductor PbTaSe 2 is in good agreement with our hrst-principles band structure calculations, 
establishing the realization of the unusual ring-shaped topological-nodal-line states in this 
compound. The topological nodal rings are protected by the reflection symmetry of the 
system. Meanwhile, the nodal rings are uniquely associated with “drumhead”-like surface 
states in a manner that resembles the connections of edge states and the nodal points in 
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graphene. Considering the nodal characteristics shared by topological nodal-line states and 
Weyl semimetals, it is very likely that the topological nodal-line materials exhibit much of 
the long-sought out exotic properties of Weyl semimetals, but with additional non-trivial 
features. For example, nodal-line states possess an extra degree of freedom for manipulat¬ 
ing novel properties of Weyl materials, which is the finite size of the nodal line. Further¬ 
more, interaction-induced instabilities that have been broadly discussed for Weyl semimetals 
should be more likely occurring in nodal-line states due to the higher density of states at the 
Fermi energy. In addition, superconductivity is induced by intercalating Pb layers to TaSe 2 , 
which also offers the essential ingredient, the Pb conducting orbitals, for the formation of 
the topological nodal-line states. Considering the intrinsic superconductivity, the spin-split 
bulk nodal-line band structure and the nontrivial surface states close to the Fermi level, it 
is possible that helical superconductivity and p-wave Cooper pairing may exist in this com¬ 
pound without the aid of the proximity effect 


25 


26| . Therefore, novel physics may arise 


from the interplay of the nodal-line states and the emergent superconductivity of PbTaSe 2 , 
which calls for further experimental investigation on PbTaSe 2 . Our ARPES measurements, 
detailed DFT simulation and theoretical analysis demonstrate the fundamental mechanism 
for realizing topological nodal-line fermions in PbTaSe 2 , and pave the way for exploring the 
exotic properties of topological nodal-line states in condensed matter systems. 


METHODS 

Single crystals of PbTaSe 2 were grown by the CVT method using chlorine in the form 
of PbCl 2 as a transport agent. For the pure synthesis of PbTaSe 2 , stoichiometric amounts 
of the elements (purity of Pb and Ta: 6N, of Se: 5N) were loaded into a quartz ampoule, 
which was evacuated, sealed and fed into a furnace (850°C) for 5 days. About 10 g of the 
prereacted PbTaSe 2 were placed together with a variable amount of PbCl 2 (purity 5N) at 
one end of another silica ampoule (length 30-35 cm, inner diameter 2.0 cm, outer diameter 
2.5 mm,). All treatments were carried out in an Argon box, with continuous purihcation 
of the Argon atmosphere resulting in an oxygen and water content of less than 1 ppm. 
Again, the ampoule was evacuated, sealed and fed into a furnace. The end of the ampoule 
containing the prereacted material was held at 850°C , while the crystals grew at the other 
end of the ampoule at a temperature of 800°C (corresponding to a temperature gradient of 





11 


2.5 K/cm) during a time of typically 1 week. Compact single crystals of sizes of up to 8 x 
5x5 mm^ were obtained. 

ARPES measurements were performed at the liquid nitrogen temperature in the beamline 
14 at the MAX-lab in Lund, Sweden. The energy and momentum resolution was better 
than 20 meV and 1% of the surface Brillouin zone (BZ) for ARPES measurements at the 
beamline 14 at the MAX-lab. Samples were cleaved in situ under a vacuum condition better 
than 1 X 10“^° torr. Samples were found to be stable and without degradation for a typical 
measurement period of 24 hours. 

STM experiments were conducted with a commercial system (Unisoku). Samples were 
cleaved at room temperature in a vucuum better than 2 x 10“^° mbar and transferred to a 
STM head precooled to 77 K. Const ant-current mode STM images were taken with chemical 
etched Pt/Ir tips. Bias voltages were applied to samples. 

We computed electronic structures using the norm-conserving pseudopotentials as im¬ 
plemented in the OpenMX package within the generalized gradient approximation (GGA) 
schemes j^, 28|. Experimental lattice constants were used j^|. A12 x 12 x 4 Monkhorst- 
Pack k-point mesh was used in the computations. The SOC effects are included self- 


consistently 


30| . For each Pb atom, three, three, three, and two optimized radial functions 


were allocated for the s, p, d, and / orbitals {s3p3d3f2), respectively, with a cutoff radius 
of 8 Bohr. For each Ta atom, d3p2d2fl was adopted with a cutoff radius of 7 Bohr. For 
each Se atom, d3p2d2fl was adopted with a cutoff radius of 7 Bohr. A regular mesh of 300 
Ry in real space was used for the numerical integrations and for the solution of the Pois¬ 
son equation. To calculate the surface electronic structures, we constructed hrst-principles 
tight-binding model Hamilton. The tight-binding model matrix elements are calculated by 


projecting onto the Wannier orbitals 


3l| . We use Pb p, Ta s and d, and Se p orbitals were 


constructed without performing the procedure for maximizing localization. 
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FIG. 1: Overview of PbTaSe2 single crystal, a, Optical image of PbTaSe 2 single crystals 
measured in this work, b, X-ray diffraction measurements showing the lattice parameters matching 
with the space group (187) P6m2. c, ARPES core level spectrum showing clear Pb 5d, Se 3d and 
Ta 4/ core level peaks, d, Resistivity as a function of temperature showing a superconducting 
transition at 3.8 K. e and f, STM topography of PbTaSe 2 indicative of a surface lattice constant 
3.2 A. 
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FIG. 2: Lattice structure and bulk bands of PbTaSe2. a, Hexagonal lattice of PbTaSe 2 . b, 
Bulk and projected (001) surface Brillouin zones, c, Calculated bulk band structure of PbTaSe 2 . 
d and e, Zoom-in band structure around H without/with the inclusion of SOC, respectively. The 
color code and arrows show the orbital component and spin orientation, respectively, f and g, 
Same as d and e but for band structure around K. 
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FIG. 3: APRES mapping and band calculation of PbTaSe2. a, ARPES spectra taken along 
M —K —r with 64 eV photons, b, DFT projected bulk bands and surface bands (bright white lines) 
of (001) surface with Pb-termination. c, ARPES isoenergy concours taken with 64 eV photons, d, 
DFT Fermi surface contour of PbTaSe 2 (001) surface. The yellow lines indicate the surface states 
on Pb-terminated (001) surface. 
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FIG. 4 : APPES measurement of PbTaSe2 and TaSe2. a, Fermi surface contour of PbTaSe 2 
(top) and TaSe 2 (bottom), b, APRES spectral cut along M — K — M of PbTaSe 2 (top) and TaSe 2 
(bottom), c, ARPES spectra of PbTaSe 2 taken with different photon energies. 
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FIG. 5: Topological-nodal rings and associated snrface states, a and b, DFT project bulk 
bands and surface bands (bright white lines) of Pb-terminated (001) surface along F —K —M and a 
generic direction K —X, respectively, c and d, Same as a and b but of Se-terminated (001) surface, 
e, The isoenergy contour showing the nodal-line states (NL) and surface states (SS). The energy is 
0.10 eV and 0.15 eV above the Fermi level for Pb- and Se- termination, respectively, as indicated 
by the red dashed lines in a-d. f and g, Bulk band structure of PbTaSe 2 with and without the 
reflection symmetry, respectively. In g, the reflection symmetry is broken by moving the Pb atom 
slightly in the vertical direction. 

















20 


FIG. 5: h, Schematic of the closed contours encircling the nodal lines for the calculation of the 
winding number. The arrows indicate the direction on the loops along which the Berry’s connection 
is integrated, j, Schematic of a spectral cut passing K. The gray shaded region indicates the 
projected bulk band associated with a single nodal ring encircling K, and the red curves depict the 
surface states (SS). 


